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Majorana fermions (MFs) are exotic particles that are their own anti-particles. Recently, the search for the 
MFs occurring as quasi-particle excitations in solid-state systems has attracted widespread interest, because 
of their fundamental importance in fundamental physics and potential applications in topological quantum 
computation based on solid-state devices. Here we study the quantum correlations between two spatially 
separate quantum dots induced by a pair of MFs emerging at the two ends of a semiconductor nanowire, in 
order to develop a new method for probing the MFs. We find that without the tunnel coupling between these 
paired MFs, quantum entanglement cannot be induced from an unentangled (i.e., product) state, but 
quantum discord is observed due to the intrinsic nonlocal correlations of the paired MFs. This finding 
reveals that quantum discord can indeed demonstrate the intrinsic non-locality of the MFs formed in the 
nanowire. Also, quantum discord can be employed to discriminate the MFs from the regular fermions. 
Furthermore, we propose an experimental setup to measure the onset of quantum discord due to the 
nonlocal correlations. Our approach provides a new, and experimentally accessible, method to study the 
Majorana bound states by probing their intrinsic non-locality signature. 

Quantum entanglement is a quantum correlation that has no classical analog. Dynamical entanglement 
has exhibited some exotic properties such as the early-stage decoherence 1 and revival 2 . Recently, it was 
found that there exists another type of nonclassical correlation coined as quantum discord 3,4 . A remark- 
able feature of the quantum discord is that it can characterize some fundamental processes of physics, including 
the quantum phase transitions 5 and the discrimination of quantum and classical Maxwell's demons 6,7 . In addi- 
tion, it can be used as a new resource for quantum information processing protocols 81 ", and be harnessed as a 
measurable quantity to implement quantum-enhanced tasks in the absence of quantum entanglement 11-13 . 

Majorana fermions (MFs) as quasi-particle excitations were anticipated to occur in different solid-state sys- 
tems, such as the 5/2 fractional quantum Hall system 14 , thep-wave superconductor 15 , and the hybrid systems of a 
topological insulator 16 . Recently, semiconductor nanowires with strong spin-orbit coupling are attracting increas- 
ing interest. In addition to the application in spin-orbit qubit 17,18 , these semiconductor nanowires provide a 
platform for demonstrating MFs 1 " 22 . Theoretically, it is shown that MFs can occur at the ends of a semiconductor 
nanowire with strong spin-orbit coupling when the nanowire is placed in proximity to an 5- wave superconduc- 
tor 19,20 . It is expected that the nanowire-superconductor system based on conventional materials is more access- 
ible to experimentally detecting the MFs 21,22 . In solid-state systems, MFs are coherent superpositions of electron 
and hole excitations. Two spatially separated MFs can form one fermionic level, which can be either occupied or 
empty, and thus defines a nonlocal qubit. The paired MFs may exhibit intrinsic nonlocal properties, so the search 
for the quantum effects induced by such nonlocal paired MFs can provide an important signature for the existence 
of MFs. Several schemes have been proposed to test the non-locality of MFs by studying the correlation behaviors 
of electron tunneling through individual MFs 23 25 . However, for the completely separated MFs, the tunneling 
events through each MF are shown to be fully independent 23 25 . Therefore, the intrinsic non-locality of the paired 
MFs cannot be demonstrated in such a context. 

A hybrid quantum system combining two or more physical systems can often utilize the strengths of different 
systems to better explore new phenomena and potentially bring about novel quantum technologies 26 . For the two 
MFs emerging at the ends of a semiconductor nanowire, because their energies are zero, it is difficult to directly 
demonstrate their intrinsic non-locality, particularly when the tunnel coupling between them becomes zero. 
Thus, to probe the emerged two MFs in the nanowire, in this work, we consider a hybrid quantum system 
consisting of two quantum dots (QDs) interacting with these two MFs. To characterize the quantum correlations 
between the two QDs mediated by the paired MFs, we use both quantum entanglement and quantum discord. We 
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Figure 1 | Schematic setup for probing quantum correlations between 
two QDs mediated by a pair of MFs. A semiconductor nanowire with 
strong spin-orbit coupling is placed in proximity to an s-wave 
superconductor, with a pair of MFs emerging at the two ends of the 
nanowire. (a) Two electrode leads are introduced to investigate the 
decoherence effects on the quantum correlations, (b) Two auxiliary QDs 
(denoted by QD L and QD K ) are introduced to measure the quantum 
correlations. The arrows indicate the tunneling directions of electrons. The 
energy levels of the QDs and the chemical potentials of the electrode leads 
are tuned by the voltages applied on them. 

show that when the paired MFs are tunnel-coupled via the wave- 
function overlapping, the maximally entangled states can be induced 
from an unentangled (i.e., product) state at specific dynamical time 
points. However, when the paired MFs are completely separated (i.e., 
no tunneling between them), quantum entanglement cannot be 
induced from this product state, but remarkably, quantum discord 
is shown to still persist owing to the intrinsic nonlocal correlation of 
MFs. This indicates that the quantum discord can indeed reflect the 
non-locality of the paired MFs. Furthermore, we propose an experi- 
mentally accessible approach to measuring the onset of quantum 
discord. Our protocol provides a new method to probe the MFs 
emerging at the two ends of the nanowire. 

Results 

Quantum dynamics. We study a hybrid structure shown in Fig. 1, 
where a semiconductor nanowire with strong spin-orbit coupling is 
placed in proximity to an s-wave superconductor. This nanowire is 
predicted to be driven into a topological superconducting phase 19,20 . 
A pair of zero-energy MFs, y 1 and y 2 , ar e anticipated to appear at the 
two ends of the nanowire. The Hamiltonian of the paired MFs is 
described by Hm = ie m y 1 y 2 /2, where e m ~e~ 1 ^ describes the tunnel 
coupling between the two MFs due to the overlap of wavefunctions. 
This inter-MF coupling damps exponentially with L (the length of 
the nanowire) and £ (the superconducting coherent length). Because 
the paired MFs have zero energies, it is difficult to directly 
demonstrate their quantum correlations, particularly when the 
tunnel coupling e m becomes zero. Thus, we introduce two QDs 
(i.e., QD] and QD 2 ) on both sides of the nanowire, each of which 
is tunnel-coupled to a MF neighboring to it. In addition, we consider 
the strong Coulomb blockade regime for each QD, so that only one 
electron is allowed therein when the tunneling process takes place. 
The Hamiltonian of the two spatially separated quantum dots is 

written as Hod = e ) d j d )' wnere d J [dj) ' s Section crea- 

tion (annihilation) operator of the dot j, with the corresponding 
energy Sj. The tunneling Hamiltonian between the QDs and the 



MFs is described by 23 " 25 H TM = ^\ _ i 2 Xj (dj — djj ty, with Xj (j = 

1, 2) denoting the tunnel coupling of the jth dot to its neighboring 
MF. 

The total Hamiltonian of the system is a sum of the above three 
parts: H sys = H M + H DD + H TM . For practical calculations, it is more 
convenient to switch from the Majorana representation to the regu- 
lar fermion one via the transformations: yi = i(f~f), and y 2 = /+/ T » 
where/is the regular fermion operator, satisfying the anticommuta- 
tion relation {f,f} = 1. After an additional local gauge transforma- 
tion di —* idi, the total Hamiltonian of the system can be written in 
the new representation as 23 25 



H sys =e m (f t /-l/2)+ J2 [¥l d i + l i{ d ]f+f d i)} 



;=l,2 



>.Ad\f+fdA+l 2 [d]fUfd 



(i) 



For simplicity, we adopt a symmetric setup with k\ = a 2 = ).. 
Moreover, the two QDs are adjusted in resonance with the two 
MFs (i.e., &\ = e 2 = 0) via the gate voltages. Also, we use ). {Ill) as 
the energy (time) unit throughout the paper. 

The QD-MF system considered in this paper includes the two QDs 
(i.e., QDi and QD 2 ) and the MFs. For the new fermion representa- 
tion, it is convenient to use the state basis \n±, n 2 , %), where n ( = 0 or 

I, with I = 1,2, M, denoting the electron occupation number in the 
left QD (i.e., QDJ, the right QD (i.e., QD 2 ) and the paired MFs, 
respectively. From Hamiltonian (1), one may notice that the electron 
numbers in both QDs and the paired MFs are not conserved, because 
a pair of electrons can be extracted out from the superconductor or 
absorbed by it. As a consequence, the parity of the isolated QD-MF 
system becomes conserved, and the Hilbert space can be split into 
two subspaces: 1 1, 1, 1), j 1, 0, 0), |0, 1, 0), and |0, 0, 1) with odd parity; 

I I, 1, 0), 1 1, 0, 1), |0, 1, 1), and |0, 0, 0) with even parity. Thus, starting 
from a product state (0, 0, 1), only odd-parity states are involved in 
the state evolution: 



|y(0) = c 1 (0|i,i,i> + c 2 (0|i,o,o) 

+ C3(0|0,l,0)-r-C 4 (t)|0,0,l) 



(2) 



With this initial condition, the time-dependent coefficients can be 
obtained, via directly solving Schrodinger equation, as follows: 

- Ae - ,n ' + Acos(Ar) - ifi sin(At) 



Ci *) = 



2A 



asin(Af) ^sin(Ar) 

^2{t)— ^ , L 3 (tJ- 



(3) 



CAt) 



Ae~ int + Acos(Ar) - ;fisin(Af) 



2A 



where A = yej^ + j 2, and Q = e m /2. From equation (2), we 

can easily obtain the density operator p(t) = | l P(t))( x I / (t)|. After 
tracing over the degrees of freedom of the MFs, the reduced density 
operator pjj) for the two QDs is obtained. This reduced density 
operator contains the nonlocal information of the paired MFs and 
serves as the starting point for calculating the quantum correlations 
(see Methods). 

In order to investigate the decoherence effect on quantum 
correlations, we also introduce two electrode leads, each weakly 
coupled to its neighboring QD. The two electrode leads can be 
treated as electron reservoirs described by the Hamiltonian 
Hieads= yi_, , e Wki c kj' where c' kj (c kj ) is the creation (anni- 
hilation) operator for the electron with wave vector k in the jth 
electrode lead and e« is the corresponding energy. The tunneling 
Hamiltonian between the electrode leads and the QDs is 
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Figure 2 | Quantum dynamics of the quantum correlations between two QDs mediated by different kinds of fermions in the absence of quantum 
decoherence. We assume a symmetric setup with Aj = X 2 = A, and Ej = £ 2 = 0- ( a ) Two QDs mediated by a pair of MFs with an inter-MF coupling 
e m = 0.5/1. (b) Two QDs mediated by a pair of MFs without an inter-MF coupling (i.e., e m = 0). (c) Two QDs mediated by a regular fermion with energy 
e c = 0. (d) Two QDs mediated by a regular fermion with e c = a, (e) Two QDs mediated by a pair of non-interacting regular fermions with energies 
e, = e 2 = 0. 



Htl= / ,j_ 12 2. /k c kj + h.c.^ . Here we employ the master 
equation approach to solving this problem. For the whole setup, 
the leads are regarded as an environment, and the central QD-MF 
device is described by a reduced density operator p(t). We assume 
that the chemical potentials of the two leads, both denoted by /.i, are 
well below the energy levels of the dots, i.e., the energy difference 8i — 
fi and e 2 — [i are much larger than the level broadening of the QDs. In 
this case, electrons can only jump outwards from the central device to 
the electrodes. In addition, we consider the regime with a weak 
coupling between the QDs and the leads. Thus, we can use the stand- 
ard Born-Markov master equation 27 . In the present case, the master 
equation can be written as 

p(t) = -i[H^,p] +r(V[d 1 ]p(t)+V[d 2 ]p(t)), (4) 

where V^p^ApA^—^^A^A^}, and the tunneling rate Y is 

associated with both the tunnel-coupling strengths and the densities 
of states of the two electrode leads. With the density operator p(t) 
solved from equation (4), we can then obtain the reduced density 
operator Pd(t) of the two QDs by tracing over the degrees of freedom 
of the MFs. 

Quantum discord versus quantum entanglement. To measure 
quantum correlations, here we use both quantum entanglement 
and quantum discord (see Methods). We show the evolutions of 
quantum entanglement and quantum discord for the two QDs 
(i.e., QDi and QD 2 ) in Fig. 2. As a typical example, we choose 
e m = 0.5l in Fig. 2a. In this case, the two MFs are tunnel-coupled, 
in addition to the intrinsic nonlocal correlation of the paired MFs. 
Equation (3) shows that there are two frequencies determining the 
oscillating behaviors of the quantum state | ^(t)). The fast frequency 
A/271 determines the period of the oscillations. From equations (2) 
and (3), it can be seen that the state of the system becomes 1,1,1) 
+ C 4 (f)|0, 0, 1) at t = nn/A (n = 1,2,3 ■ ■ ■ ), i.e., the state of the two 
QDs is just an entangled state Ci(f)|l, 1) + C 4 (0|0, 0). As shown in 



Fig. 2a, the quantum entanglement peaks at t = nn/A. The slow 
frequency Q/2n modulates the height of oscillating peaks and 
determines the period of the modulating cycles. It can be seen 
that the quantum discord exhibits similar behaviors of oscillations 
except for the very short initial and final periods of each cycle, where 
the quantum entanglement becomes close to zero. Also, as shown in 
Fig. 2a, the quantum entanglement (discord) can reach its maximal 
value Cd = 1 (A; = 1) at some specific time points. When e m = 0, 
which corresponds to the long-length limit of the nanowire, the two 
MFs are completely separated in space. In this case, there is no direct 
tunnel coupling between the two MFs, but their intrinsic nonlocal 
correlation still exists. As shown in Fig. 2b, the temporal behaviors 
of quantum entanglement and quantum discord are very different at 
e m — 0. While the quantum entanglement remains zero with time, 
the quantum discord oscillates with a smaller amplitude. Previous 
studies have found it hard to demonstrate the intrinsic nonlocal 
correlation between the two MFs in the electrical measurement 
signals, including the electrical current and the current noise 23 25 . 
Thus, it will be important to select a physical quantity that can 
significantly reflect the non-locality of the MFs. Interestingly, our 
results reveal that the quantum discord can clearly serve as an 
indicator showing the intrinsic nonlocal correlation of the MFs in 
contrast to the insensitive quantum entanglement considered in our 
system. 

For the purpose of discriminating MFs from regular fermions, 
we also show the quantum correlations of the two QDs (i.e., QDi 
and QD 2 ) induced by regular fermions that couple to them. We 
first consider the case with the two QDs coupled by a common 
regular fermion. Different from equation (1), the Hamiltonian 
of this QD-fermion system is given by H sys = e c c'c+ 

Ejdjdj + Aj (djc+c^dj\ J , where c 1 (c) is the electron creation (anni- 
hilation) operator of the regular fermion with energy e c , and kj (j = 1, 
2) denotes the tunnel coupling between the jth dot and the regular 
fermion. We assume that the two QDs is initially empty, but the 
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fermionic level e c is occupied by a single regular fermion such as the 
fermionic quasi-particle excitation in the nanowire. When the ferm- 
ionic level e c is resonant to the energy levels e t of the two QDs, both 
quantum discord and quantum entanglement of the two QDs 
induced by the fermion have identical dynamical evolutions (see 
Fig. 2c); when e c is off-resonant to £ ; , these quantum discord and 
quantum entanglement have different amplitudes, but still exhibit 
very similar temporal behaviors. Moreover, we consider the case 
involving a pair of non-interacting regular fermions, where each 
fermion couples to a QD. The Hamiltonian of this QD-fermion 

system is H sys = , \ejdj dj + EjC- Cj + lj(d'j Cj + cj dj\\> where 

cj (cj) is the electron creation (annihilation) operator of the regular 
fermion with energy ej, and Xj (j = 1, 2) denotes the tunnel coupling 
between the jth quantum dot and the regular fermion that couples to 
it. Also, we assume that the two QDs are initially empty, but each 
fermionic level e t is occupied by a single regular fermion. Because 
there is no intrinsic non-locality between these two non-interacting 
regular fermions, both quantum discord and quantum entanglement 
of the two QDs remains zero with time (see Fig. 2e). As compared 
with Figs. 2a and 2b, these very different temporal behaviors in 
Figs. 2c-2d can be used to discriminate the MFs from the regular 
fermions. 

We also study the temporal behaviors of both quantum entangle- 
ment and quantum discord for the QD-MF system in the presence of 
two electrode leads. Here we choose a tunneling rate Y = 0.051 to 
ensure a weak coupling between the central QD-MF system and the 
two electrode leads. These two leads play the role of a fermionic 
environment which induces decoherence on the evolutions of 
quantum correlations. With tunnel coupling between the paired 
MFs, e.g., e m =Q.5l in Figs. 3a, both quantum entanglement and 
quantum discord can be generated from the product state, but the 
oscillating quantum entanglement decays faster with time than the 
quantum discord. This implies that the quantum discord, as a new 
resource of quantum information processing, is more robust than the 
quantum entanglement in this model system. Moreover, without 
tunnel coupling between the paired MFs, i.e., e m = 0, the quantum 
entanglement still remains zero with time, but the oscillating 
quantum discord decays due to the decoherence induced by the 
electric leads (see Fig. 3b). 

A measurement scheme. Note that equation (2) can be rewritten as 
|«P(t)) = b 1 (t) e *W|l,l,l) + fe(t)|l,0,0) 

(5) 

+ fo 2 (0|0,l,0) + fc 3 (0e* () |0,0,l), 

where the real parameters £>i(f)> b 2 (t), b 3 (t), (j>i(t) and <f> 3 (t) are 
determined by C,(f)> i = 1, 2, 3, and 4, in equation (3). After 
tracing over the degrees of freedom of the MFs, the reduced 
density operator of the two QDs (i.e., QDj and QD 2 ) can be 
obtained via equation (5) as 



with A(j){t) = ^(r) - <j) 3 (t). Here b\(t), b\(t) and b 2 3 (t), i.e., the 
diagonal matrix elements in equation (6), correspond to the 
probabilities of both dots occupied, either dot occupied, and both 
dots empty, respectively. These probabilities can be directly obtained 
by joint measurements on the electron occupation of QD : and QD 2 
(see, e.g., Ref. 30). In order to obtain the values of the phase difference 
A(/)(t), we introduce, as Fig. lb shows, two auxiliary QDs (denoted by 
QD L and QD R ) adjacent to QD t and QD 2 , respectively. First, instead 
of |0), these two auxiliary QDs are initially prepared in the 



superposition state: yji 



L(R) 



= (|0) z 





0 


bi[t)b 3 [t)(P*®\ 




0 b\{t) 


m 


0 


,(6) 


o bl(t) 


m 


0 




0 


m ) 





/L(8)T|i/i (8 ) j j Vi Second, the 
QD L (QD K ) is adjusted in resonance with the QDj (QD 2 ) (i.e., e L 
= 8i = 0, and e R = e 2 = 0) at the moment r, via tuning the gate 
voltages. Simultaneously, the tunnel coupling between QD l (QD 2 ) 
and its adjoining MF is switched off by tuning the voltages of the 
electrical gates that control the tunnel barrier between them. The 

Hamiltonian of the total system becomes Htd = T\ (d\di + d\ dij + 

T 2 (d' R d 2 + dldg^j , where T 1 (T 2 ) is the tunnel coupling between QD L 

(QDj;) and QD] (QD 2 ) and assume T l = T 2 = T. Therefore, the 
system further evolves under the Hamiltonian H TD , from the initial 
product state | ip L ) ® | ip R ) ® ^(f)). Third, along with the evolution 
of the total system, we choose a specific moment, t + At, to jointly 
measure the occupation probabilities of both QD : and QD 2 , i.e., P n , 
Pio> Poi> an d Poo which denote the probabilities of both dots 
occupied, single QD t occupied, single QD 2 occupied, and both 
dots empty, respectively. When At = iz/iT, these probabilities can 
be obtained as 



P u W = - [9bl(t) + Sb l 2 (t) + bl(t)+2b 1 (t)b 3 (t) cos Aftt)] , 
Poo(0 = y 6 [b 2 l (t)+8b 2 2 (t)+9b 2 3 (t) +2^(^(1) cos A(j>[t)] , 
V 10 (t) = ~[3b 2 (t) + 8b 2 2 (t) + 3b 2 3 (t)-2b 1 (t)b 3 (t) cosA^(f)], 



(7) 



and P O i(0 = PioW- This measurable occupation probabilities 
include the phase difference A0(r), in addition to b^t), b 2 (i), and 
b 3 {i) that have previously been obtained via joint measurements as 
well. Thus, it is experimentally feasible to obtain the reduced density 
operator p d (t) of the two QDs mediated by the paired MFs. Once the 
reduced density operator pa(t) is determined experimentally, one can 
then derive the quantum entanglement and the quantum discord of 
the two QDs using this reduced density operator (see Methods). 




10 

Time (VX) 

Figure 3 | Quantum dynamics of the quantum correlations between two 
QDs mediated by a pair of MFs in the presence of quantum decoherence. 

We assume a symmetric setup with Aj = ~a 2 = I, and £\ = £ 2 = 0. The 
tunneling rate between the central QD-MF device and the two electrode 
leads is chosen as T = 0.05A. (a) e m = 0.5k, and (b) e m =0. 
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Discussion 

We have investigated the dynamics of quantum correlations between 
two QDs mediated by a pair of MFs. We find that in the presence of 
tunnel coupling between the two MFs, maximally entangled states of 
the two QDs can be generated from an unentangled state under 
proper dynamical conditions. Notably, in the absence of tunnel 
coupling between MFs, while the quantum entanglement remains 
zero, quantum discord can be induced from this unentangled state 
owing to the intrinsic nonlocal correlations of the MFs. While the 
conventional methods fail to reveal the non-locality information of 
the MFs, our demonstrated new feature indicates that the quantum 
discord can serves as a novel quantum signature reflecting the 
intrinsic non-locality of the MFs. Physically, it shows that quantum 
discord captures some interesting correlation features of the system 
that is insensitive to entanglement. Also, we show that these features 
do not exist when replacing MFs with regular fermions. Therefore, 
they can be used to discriminate the MFs from the regular fermions. 
Furthermore, we have proposed a scheme for experimentally deter- 
mining the quantum discord of the two QDs mediated by the paired 
MFs. In addition, we further investigate the decoherence effects on 
the quantum correlations, and the results show that the quantum 
discord is more robust than the quantum entanglement. Because the 
quantum discord is a measurable physical quantity and the intrinsic 
nonlocal nature is an extraordinary character of the MFs, our 
approach provides a new, and experimentally feasible, method to 
probe the MFs by demonstrating their intrinsic non-locality via the 
quantum discord. 

Methods 

With the reduced density operator p d of the two QDs (i.e., QD : and QD 2 ), the 
concurrence 28 , as a measure of the quantum entanglement, can be obtained as 

Qj= maxjcv^-^V!- \f^ — v^}, (8) 

where A 2 , A 3 and A 4 are the square roots of the eigenvalues in decreasing order for 
the matrix 

G=P d (0 y ®0 y )p* d {0y®0y), (9) 

where a y is a Pauli matrix, and p* d is the complex conjugate of p d . 

Quantum discord is a measure of quantum correlations obtained by subtracting 
the classical correlations from the total correlations 3,4 . For a bipartite quantum 
system including subsystems A and B, the total correlations are defined, in terms 
of quantum mutual information, as I(p A : p B ) = S(p A ) + S(p B ) — S(p AB ), where 
p AB denotes the density operator of the bipartite system, with p A (p B ) being the 
reduced density operator of the subsystem A (B), and S(p AB ) = — Tr(p AB log 2 PAn) 
is the entropy of the bipartite quantum system. The total classical correlations 
based on the POVM measurement of the subsystem B are defined as 3,29 
/(PabK- 8 *}) = s (p A ) - mia {B t } ^2 k PkS(p k ), where {B k } is a complete set of pro- 
jectors of local measurements performed on the subsystem B\ p k = Tr[(7 ® 
B^PabO ® Bk)] denotes the probability of obtaining the measurement outcome k, 
with 1 being the identity operator for the subsystem A; p k ~ (I ® B^)p AB {l (x) B^)j 
pk is the state of the bipartite system, conditioned on the measurement outcome 
labeled by k. The quantum discord is defined as 3,4 

D B (p AB ) = min tBl) {l( P A : /) -J(p AB \{B t })] . (10) 

For our system, the A and B parts are the two QDs {i.e., QDi and QD 2 ), and p AB 
— pd- To compute classical correlations, we adopt the specific case of von 
Neumann local measurements, which are provided by the orthogonal projectors 
Bi c — RYli c R li , where {11^} — {j/c){/cj} defines the computational basis R, with k = 

0, 1 corresponding to the empty and occupied states, respectively. As a rotation 
operator, V can be parameterized as 

/ cosf sinfe-'A 
R= [-sml^ -cosf j (U) 

where 0 £ [0, n], and (j) E [0, 2n\. Hence, the quantum discord given in equation 
(10) can be directly obtained by minimizing over all angles 9 and </>. 
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